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Several key problems of QCD sum rules in the spin-0 glueball channels are 
resolved by implementing nonperturbative short-distance physics from direct 
instantons and topological charge screening. A lattice-based instanton size dis- 
tribution and the IR renormalization of the nonperturbative Wilson coefficients 
are also introduced. Results of a comprehensive quantitative sum rule analysis 
are reviewed and their implications discussed. 



1. Introduction 

The glucinium states of QCD have remained intriguing for almost four 
decades 1 . Their "exotic" nature reflects itself not least in several long- 
standing problems which the QCD sum rule approach faces in the spin- 
glueball channels 2 . In the scalar (0 ++ ) glueball correlator, in particular, 
the departure from asymptotic freedom sets in at unusually small distances 3 
and the perturbative Wilson coefficients of the standard operator product 
expansion (OPE) proved inadequate to establish consistency both among 
the ++ glueball sum rules and with an underlying low-energy theorem 4 . Al- 
though nonperturbative contributions due to direct (i.e. small) instantons 5 
were early candidates for the missing short-distance physics 6 , insufficient 
knowledge of the instanton size distribution 5 ' 7 prevented their quantitative 
implementation at the time. Only recently, the derivation of the exact in- 
stanton contributions (to leading order in h), their duality continua and 
the corresponding Borel sum-rule analysis 4 showed that direct instantons 
indeed solve the mentioned key problems in the scalar glueball channel. Be- 
low we will outline a more thorough and systematic treatment 9 which elimi- 
nates artefacts of earlier approximations, significantly modifies the sum-rule 
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results in both spin-0 channels and improves their reliability. Moreover, we 
implement for the first time topological charge screening contributions to 
the OPE in the h channel 9 and show how those resolve earlier sum- 
rule inconsistencies 10 . Finally, we review several new predictions for spin-0 
glueball properties. 



2. Correlators and sum rules 

The spin-0 glueball sum rules are based on the scalar (0 ++ ) and pseu- 
doscalar (0~ + ) correlations functions 

Il G (-q 2 )=i J d 4 xe^ (0\TO G (x)O G (0)\0) (1) 

where Oq with G e {S, P} are the gluonic interpolating fields (of low- 
est mass dimension) O s = a s G a ^G a ^ and O p = a s G£„G a ^ with 
Gpv = (1/2) e^upaGpa (the definition in Ref. 9 contains a typo). The zero- 
momentum limit of the correlator (1) is governed by the low-energy theo- 
rems (LETs) 

n s (<z 2 =0) = ^<aG 2 ) (2) 
Oo 

in the scalar 11 and (for three light flavors and m u ^ <C m s ) 

U P ^=0)={8n) 2 ^^-(qq) (3) 
m u + m d 

in the pseudoscalar channel 12 . (Note that Eq. (3) vanishes in the chiral 
limit.) Consistency with the low-energy theorems places stringent con- 
straints on the sum rules which cannot be satisfied without nonperturbative 
short-distance physics 4 ' 9 . 

Contact with the hadronic information in the glueball correlators is 
established by means of the dispersive representation 

where the necessary number of subtractions is implied but not written ex- 
plicitly. The standard sum-rule description of the spectral functions 

Im n<? h) (s) = Im n^°' e) (s) + Im n^ mt) (a) (5) 

contains one or two resonance poles in zero-width approximation, 
l m Y[iP° le ) ( s ) = irJ2i = i fGi m Gi& ( s ~ m Gi)> anal tnc local-duality contin- 
uum ImllQ™'' (s) = 8 (s — so) Imlfg ^' (s) from the IOPE discontinu- 
ities in the "duality range" which starts at an effective threshold sq. 
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In order to write down QCD sum rules, the Borel-transformed IOPE 
and dispersive representations of the correlators - weighted by k powers of 
— Q 2 - are matched in the fiducial r-region (r is the Borel variable) and 
rearranged in terms of the continuum-subtracted Borel moments 1Za,k as 

TZ G . k (r; a ) = - dss k Imng OPB) (a) e~ ST (6) 
h 

i=i 

The pole contributions of interest (and for k = — 1 the crucial subtraction 
terms) are then isolated on the RHS, and the hadronic parameters m Gi , f Gi 
and So can be determined numerically. 



3. IOPE 

Our theoretical framework, the instanton-improved operator product ex- 
pansion (IOPE), factorizes the correlators at large, spacelike momenta 
Q 2 = —q 2 3> A^qcd into contributions from "hard" field modes (with 
momenta |fc| > p) in the Wilson coefficients Cd (Q 2 ) an d "soft" field 
modes (with \k\ < ji) in the "condensates" (Od^ of operators Ojj with 

increasing dimension D. Previous glucball sum rules based on the OPE 
with purely perturbative Wilson coefficients were plagued by notorious in- 
consistencies between the predictions of different moment sum rules and 
by massive LET violations. Moreover, the soft nonperturbative condensate 
contributions were exceptionally small. 

We have therefore analyzed hard nonperturbative contributions to the 
Wilson coefficients. They are strongly channel dependent and due to direct 
instantons and topological charge screening. The instanton contributions 
to the spin-0 coefficients 9 , 

4 /+J) (* 2 ) ^jdpn (p) 1 2Fl ( 4 , 6, §,-£), (7) 

are large (while those in the 2 ++ tensor channel vanish) and add to the unit- 
operator coefficients Cq . The imginary part of their Fourier transform at 
timclike momenta generates the Borel moments 4 

Tl{ 1+I) (t) = -2 7 7r 2 ( 5 fc ,_ 1 / dpn (p) - 2V f dp (8) 



x n (p) p 4 / dss k+2 J 2 (sfsp) Y 2 (Tip) e 
Jo 
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which are similar or larger in size than the perturbative ones. The evaluation 
of these moments requires as the sole input the (anti-) instanton distribu- 
tion n (p) which is implemented by means of a lattice-based Gaussian-tail 
parametrization with the correct small-/9 behavior, 



> 18 n ( p\ 4 ( 2 6 p 



(for N c = Nf = 3), which was introduced in Rcf. 9 and shown to pre- 
vent several artefacts of the oversimplified "spike" approximation n(p) = 
nS (p — p) on which all previous direct-instanton calculations had relied. 

As an additional benefit, the realistic size distribution allows for a gauge- 
invariant IR renormalization which excludes large instantons with size p > 
p^ 1 from the Wilson coefficients, 

n{p) -^h(p;p) = e (p- p,- 1 )n(p) (10) 

(8/3 is a soft step function). The instanton-induced p dependence turns out 
to be relatively weak for p < as necessary to compensate its per- 

turbative counterpart. Neglect of this renormalization, although common 
practice in perturbative Wilson coefficients, would significantly contaminate 
the results, e.g. by missing the reduction of the direct-instanton density 

/>oo />oo 

h= dpn(p) — > / dpn(p: p) = n(p) . (11) 
Jo Jo 

Another important renormalization effect is the reduction of the instanton 
contributions to the pseudoscalar relative to the scalar sum rules. 

The instanton's self-duality causes a strongly repulsive contribution to 
the h channel, with seemingly detrimental impact on the sum rules 10 : 
the glueball signal disappears and both unitarity and the LET (3) are badly 
violated. The origin of these problems can be traced to the neglect of topo- 
logical charge screening 9 . Due to their high channel selectivity and a small 
screening length Ad ~ m^, 1 ~ 0.2 fm, the model-independent screening 
correlations 14 

n^r ' (x) ~ -2 V f ( m (x) m (0)) (12) 

(£ is the overall topological charge density (= n for instantons) and 770 the 
flavor-singlet part of the 77') affect almost exclusively the h Wilson co- 
efficients. They arise from the axial-anomaly induced attractive (repulsive) 
interaction between topological charge lumps of opposite (equal) sign due 
to 770 exchange 13 and, after correcting for 770 — r/s mixing, add the terms 

< r) W = (|r + + F^e-<- + F*mfe-<T (13 ) 
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to the pseudoscalar IOPE moments. 

They strongly reduce the direct-instanton induced repulsion and resolve 
the disastrous problems mentioned above 9 : positivity of the spectral func- 
tion is restored, the four h Borcl sum rules (6) are stable and contain 
consistent pseudoscalar glueball information. (All previous analyses had 
discarded the k = — 1 sum rule and thereby missed valuable first-principle 
information and LET consistency checks.) 

4. Results and discussion 

We have implemented direct instanton and topological charge screening 
contributions into the OPE coefficients of the spin-0 glueball correlators, 
evaluated their duality continuua and demonstrated how these contribu- 
tions resolve the problems of previous QCD sum rule analyses 9 . A lattice- 
based instanton size distribution and the gauge-invariant 1R renormaliza- 
tion of the nonperturbative Wilson coefficients were also introduced. Quark 
admixtures, and thereby quarkonium mixing effects, enter through quark 
loops, the instanton size distribution and the condensates. 

In the scalar channel, the sizeable direct instanton contributions are in- 
dispensable for mutually and LET consistent sum rules. Their improved 
treatment reduces our earlier (spike-distribution based) result for the ++ 
glueball mass to rag = f .25 ± 0.2 GeV. (The mass stays well beyond f 
GeV, however, in contrast to obsolete predictions based on purely perturba- 
tive coefficients.) This value is somewhat smaller than the quenched lattice 
results 15 (which will probably be reduced by light-quark effects) and con- 
sistent with the broad glueball state found in a comprehensive if -matrix 
analysis 16 . The systematics among our different Borel moments likewise 
indicates a rather large width of the scalar glueball, Lg > 0.3 GeV. 

Our prediction for the scalar glueball decay constant, fs = 1 .05 ± 0.1 
GeV, is several times larger than the value obtained when ignoring the non- 
perturbative Wilson coefficients. This implies an exceptionally small ++ 
glueball size, in agreement with several lattice results 17 . Another stringent, 
OPE- and sum-rule-indcpcndcnt consistency check of the instanton con- 
tributions and their fs enhancement provide numerical simulations in an 
instanton ensemble 18 (ILM). For \x\ < 0.5 fm the scalar ILM and IOPE 
correlators are indeed very similar, and the ILM prediction f^ LM ^ = q g 
GeV is similarly large. This indicates a robust instanton effect and rules 
out that the large fs "may signal some eventual internal inconsistencies 
in the treatment of the instanton contributions" 19 . A subsequent model- 
independent confirmation of the fs enhancement was supplied by the first 
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direct (quenched) lattice result 20 fg — 0.86 ± 0.18 GeV which is consis- 
tent both with our prediction and the ILM value (the latter is practically 
independent of quark quenching 18 ). 

Our prediction for fg implies enhanced partial widths for radiative J/ ip 
and T decays into scalar glueballs and is therefore relevant for experimental 
glueball searches, e.g. in the CLEO and BES data on T — > 7/0 and other 
decay branches. Since the exceptionally small size and large decay constant 
are particular to the scalar glueball and since only part of the scalar decay 
constant contributes to the radiative production rates 21 , however, the above 
results are not ruled out by experimental data on decays into h glueballs 
(cf. Ref. 19). 

In the pseudoscalar (0 '") glueball correlator we have identified and im- 
plemented a new type of nonperturbative contributions to the Wilson coef- 
ficients, due to topological charge screening. Roughly speaking, the screen- 
ing effects "unquench" the direct instanton contributions, thereby restoring 
unitarity, the axial Ward identity and the resonance signals. Consistency 
among all moment sum rules and with the underlying LET is also achieved, 
and the resulting mass prediction rap = 2.2 ±0.2 GeV lies inside the range 
of quenched and unquenched lattice data. The coupling fp = 0.6 ± 0.25 
GeV is somewhat enhanced by the topological short-distance physics, af- 
fecting radiative production rates, the 77 — > Gpir cross section at high 
momentum transfers and other glueball signatures. 

The crucial impact of the nonperturbative Wilson coefficients on both 
spin-0 glueball correlators is particularly evident in the interplay between 
their subtraction constants. Indeed, the notorious consistency problems 
which plagued previous ++ glueball sum rules were primarily caused by the 
large LET-induced subtraction constant Tig (0) ~ 0.6 GeV 4 (cf. Eq. (2)): it 
cannot be matched by perturbative Wilson coefficients (since n^ er ^ (0) = 

0) and requires the direct instanton contribution dgjp (0) = ±2 7 TT 2 n ( ur — 
0.63 GeV 4 . At first sight this seems to imply a conflict with the much 
smaller LET subtraction constant lip (0) ~ -0.02 GeV 4 (cf. Eq. (3)) in 
the H channel, however, since the instanton contributions to both spin-0 
correlators are equal (up to a sign). Here the topological screening contribu- 
tions ITp Cr) (0) ~ 0.59 GeV 4 from Eq. (13) prove indispensable: they restore 
consistency by canceling most (and in the chiral limit, where lip (0) — > 0, 

all) of the instanton contributions: U { p Cr) (0) + Up 1 ' 1 ' (0) ~ lip (0). 

To summarize: contrary to naive expectation, the nonperturbative con- 
tributions to the OPE of the spin-0 glueball correlators reside primarily in 
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the Wilson coefficients (i.e not in the condensates) and are closely related 
to the topological vacuum structure. The nonperturbative short-distance 
physics resolves long-standing consistency problems of the associated QCD 
sum rules and generates a rather diverse set of new glueball predictions. A 
large part of the ++ glueball mass and binding originates from direct in- 
stantons, for example, while their net effects in the CP channel are smaller 
and more subtle, due to cancellations between instanton and topological 
charge screening contributions. 

This work was supported by FAPESP and CNPq of Brazil. 
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